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Abstract 

It has been long conjectured that the crossing numbers of the complete bipartite 
graph K m ^ n and of the complete graph K„ equal Z{m,n) := |_fj L 22 ^"] LtJ \ 
and Z{n) := \ [§J \ ^^\ > respectively. In a 2-page drawing of a graph, the 

vertices are drawn on a straight line (the spine), and each edge is contained in one of the 
half-planes of the spine. The 2-page crossing number v 2 (G) of a graph G is the minimum 
number of crossings in a 2-page drawing of G. Somewhat surprisingly, there are 2-page 
drawings of K m , n (respectively, K n ) with exactly Z(m, n) (respectively, Z{n)) crossings, 

thus yielding the conjectures (I) v 2 (K m ^ n ) = Z(m,n) and (II) u 2 {K n ) = Z(n). It is 
known that (I) holds for min{m,n} < 6, and that (II) holds for n < 14. In this paper 
we prove that (I) holds asymptotically (that is, \m\ rL ^ f00 V2{K myn )/Z(m ) n) — 1) for 
m = 7 and 8. We also prove (II) for 15 < n < 18 and n 6 {20, 24}, and establish the 
asymptotic estimate 

lim u 2 (K n )/Z(n) > 0.9253. 
The previous best-known lower bound involved the constant 0.8594. 
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maximum cut, Goemans- Williamson max-cut bound 
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1 Introduction 

We recall that the crossing number cr(G) of a graph G is the minimum number of pairwise 
intersections of edges (at a point other than a vertex) in a drawing of G in the plane. Besides 
their natural interest in topological graph theory, crossing number problems are of interest 
because of their applications, most notably in VLSI design [25] , 
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Also motivated by applications to VLSI design, Chung, Leighton and Rosenberg [5] studied 
embeddings of graphs in books: the vertices are placed along a line (the spine) and the edges 
are placed in the pages of the book. In a book drawing (equivalently, k-page drawing, if the 
book has k pages), crossings among edges are allowed. The k-page crossing number Uk{G) 
of a graph G is the minimum number of crossings of edges in a /c-page drawing of G. 

Clearly, a graph G has V\(G) = if and only if it is outerplanar. Closely related to 1-page 
drawings are circular drawings, in which the vertices are placed on a circle and all edges 
are drawn in its interior. It is easy to see the one-to-one correspondence between 1-page 
drawings and circular drawings. 

In a similar vein, 2-page drawings can be alternatively modelled by drawing the vertices of 
the graph on a circle, and imposing the condition that every edge lies either in the interior 
or in the exterior of the circle (see Figure [T|). In this paper we shall often use this equivalent 
circular model for 2-page drawings, as well as the usual spine model. It is known that the 
family of graphs G with i>2{G) = is precisely the family of subgraphs of Hamiltonian planar 
graphs [3]. As a consequence, there exist planar graphs G with ^(G) > 0, in contrast to 
the case of the normal crossing number. In fact, it was shown that all planar graphs may 
be embedded without crossings in 4-page books, and that four pages are necessary [3D] . 



Figure 1: A 2-page drawing of (a) in the spine model; and (b) in the circular model. 

Masuda et al. [281 [29] proved that the decision problems for v\ and U2 & re NP-complete. 
Shahrokhi et al. [34j gave an approximation algorithm for u^(G), as well as applications to 
the rectilinear crossing number. A more recent, additional motivation for studying A;-page 
crossing numbers comes from Quantum Dot Cellular Automata [36] . 

Several interesting algorithms and heuristics have been proposed for producing 1- and 2- 
page drawings (see for instance [U [181 Q21 1201 EI])- As with the usual crossing number, 
the exact computation of Vk{G) (for any integer k) is a very challenging problem, even 
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for restricted families of graphs. In this direction, Fulek, He, Sykora, and Vrt'o [9], He, 
Salagean, and Makinen [17], and Riskin [33J have computed the exact 1-page and 2-page 
crossing numbers of several interesting families of graphs. 

1.1 Drawings of K m ^ n and K n 

Turan asked in the 1940's: what is the crossing number of the complete bipartite graph 
K m>n ? There is a natural drawing of K m ^ n with exactly Z(m,n) := |_|J \_ J1: ^-\ [_?rj L^TT"] 
crossings (see Figure [2]), and so cr(K m)n ) < Z(m,n). 




Figure 2: A drawing of K$ t Q with Z(5, 6) = 24 crossings. By performing a homeomorphism 
from the plane to itself that takes the dotted curve to a straight line, the result is a 2-page 
drawing of K^^ with the same number of crossings. 

Perhaps the foremost open crossing number problem is Zarankiewicz's Conjecture, dating 
back to the early 1950's [H]: 

cr(if min ) = Z{m,n). (1) 

This conjecture has been verified only for min{m,n} < 6 [22], and for the special cases 
(m, n) G {(7, 7), (7, 8), (7, 9), (7, 10), (8, 8), (8, 9), (8, 10)} M- 

On a parallel front, there are drawings of the complete graph K n with exactly Z(n) := 
I L§J l_^T^"J L^ir^J l_^^J crossm g s (f° r every n), and so cx{K n ) < Z(n). These drawings 
inspired the still open, long-standing Harary-Hill Conjecture |15| : 

cr(^ n ) = Z(n). (2) 
This conjecture has been verified for n < 12 [30J. 

For a detailed account on the history of ([T]) and ([2]), we refer the reader to the lively survey 
by Beineke and Wilson [2]. 
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1.2 2-page drawings of K m ^ n and K n 

The drawing in Figure [2] is easily generalized to yield a drawing of K m ^ n with Z(m, n) 
crossings. As mentioned in the caption of this figure, such a drawing is easily transformed 
into a 2-page drawing of K m ^ n with the same number of crossings. Thus, there exist 2-page 
drawings of K m ^ n with Z(m, n) crossings. 

On the other hand, it is somewhat surprising that there exist 2-page drawings of K n with 
exactly Z{n) crossings, for every positive integer n ([H]; see also [16]). 

These observations imply that vy,{K m ,n) < Z(m,n) and U2(K n ) < Z(n). Since obviously 
cr(G) < 2*2 (G) for every graph G, ([1]) and ([2]) immediately imply the following conjectures: 

V 

V2{K m>n ) = Z(m,n). (3) 
u 2 (K n ) = Z{n). (4) 

Even though ([3|) and (j4|) are (at least in principle) weaker than the corresponding ([I]) and 
([2]), and even though the 2-page crossing number problem can be naturally formulated in 
purely combinatorial terms, our current knowledge (prior to this paper) on (J3|) and (J4|) is 
not substantially better than our knowledge on ([1]) and ([2]). Indeed, the only step ahead is 
the proof by Buchheim and Zheng [4] that 1^2(^13) = Z(13) (from which a routine counting 
argument yields that ^2(^14) = Z(14)). The best general lower bounds known for V2{K mjTl ) 
and V2{K n ) are the same as those known for cr(K mtn ) and cr(iT n ), and the same is true for 
the asymptotic ratio lim n _ s>00 V2(K n )/Z(n), whose best current estimate is exactly the same 
as the asymptotic ratio lim n _^ 00 cr(K n )/Z(n), namely 0.859 [23J. 



1.3 Main results 

Our main results in this paper offer a substantial improvement on our knowledge of V2{K m ,n) 
and V2{K n ) over our knowledge of ci{K m , n ) and cr(K n ). 

Theorem 1. The 2-page Harary-Hill Conjecture holds for all m < 18 and for m = 20 and 
24: 

v 2 {K m ) = Z(m) for all m < 18 and for mG{20,24}. (A) 

Moreover, the asymptotic ratio between the 2-page crossing number of K n and its conjec- 
tured value satisfies: 

Jim V -2^A > 0.9253. (B) 

n->oo Z(n) 

Theorem 2. The 2-page Zarankiewicz's Conjecture holds in the asymptotically relevant 
term for m = 7 and 8. That is: 

V2(K 7 , n ) = (9/4)n 2 + O(n) = Z(7,n) + 0{n), and 

^ 2 (^8,n) = 3n 2 + 0(n) = Z(8, n) + 0(n). 
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Therefore, 

\hn°4^4 = l and lim = 1. 

n->oc Z(7,n) n->oo Z(8,n) 

Outline of this paper 

The rest of this paper is structured as follows. In Section [21 we review the reformulation 
(first unveiled by Buchheim and Zheng [3]) in which the problem of calculating V2{K n ) is 
shown to be equivalent to a maximum cut problem on an associated graph G n . In Section 3 
we invoke a result by Goemans and Williamson that provides an upper bound on the size of 
the maximum cut of a graph; this bound may be computed via semidefinite programming. 
Using these ingredients, in Section @] we present the numerical computations that establish 
Theorem [TJ In Section [5] we formulate a quadratic program whose solution yields a lower 
bound on V2{K m ,n)- In Section [6] we analyze the semidefinite programming relaxation of 
this quadratic program, and in Section [7] we give the numerical computations that prove 
Theorem [2j In Section [8] we present some concluding remarks. 

2 Formulating v^iK-n) as a maximum cut problem 

Buchheim and Zheng [4] unveiled a natural reformulation of the fixed linear crossing number 
problem (FLCNP) as a maximum cut problem. Their results imply, in particular, that 
V2{K n ) can be obtained by computing the maximum cut size in a certain graph G n — 
(V n ,E n ), with V n and E n defined as follows. 

Consider a Hamiltonian cycle with vertices V\,V2, ■ ■ ■ ,v n . Let V n be the set of chords of the 
cycle, that is, the edges ViVj with Vi and Vj at cyclic distance at least 2. Now to define E n , 
let two chords v^Vj and VkVe be adjacent if they intersect. This construction is illustrated 
in Figure [3] for n = 5. 




Figure 3: The chords v\v?, and V2V5 form adjacent vertices in the graph G5. 
Thus \V n \ = (2) - n, and it is easy to check that \E n \ = (^). The automorphism group of 
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G n is isomorphic to the dihedral group D n , and there are d — 1 orbits of vertices, where 
d = \n/2\ . The equivalency classes of vertices (i.e. vertices belonging to the same orbit) may 
be described as follows: since vertices correspond to chords in P n , the chords that connect 
vertices of P n at the same cyclic distance belong to the same equivalency class. The vertices 
corresponding to chords with cyclic distance i have valency i(i — 1) + 2(i — l)(d — i), as is 
easy to check. 

Now for a graph G = (V, E) and a subset W C V, cutn/(G) denotes the number of edges 
with precisely one endpoint in W, and maxcut(G) is the maximum value of cutw(G) taken 
over all subsets W C V. 

The next lemma follows immediately from Theorem 1 in [4]. We sketch the proof for the 
sake of completeness. 

Lemma 3. 

vz{K n ) = \E n \- maxcut(G n ). 

Proof. Given a two page (circle) drawing of K n , define W C V n as the chords that are 
drawn inside the circle. The edges of E n with precisely one endpoint in W now correspond 
to edges of K n that do not cross in the drawing. □ 

As a consequence of this lemma, one may calculate v^iKn) for fixed (in practice, sufficiently 
small) values of n by solving a maximum cut problem. This was done by Buchheim and 
Zheng [1] for n < 13, by solving the maximum cut problem with a branch-and-bound 
algorithm (Bucheim and Zheng applied the technique to many other graphs as well). Using 
the BiqMac solver |31| . we have computed the exact value of viiKn) for n < 18 and for 
n G {20,24} (statement (A) in Theorem [TJ see Section U]). 



3 The Goemans-Williamson max-cut bound 



We follow the standard practice to use W xq (respectively, C pxq ) to denote the space of 
p x q matrices over R (respectively, C). For A £ R pxp , the notation A y means that 
A is symmetric positive semidefinite, whereas for A S C pxp , it means that A is Hermitian 
positive semidefinite. 

Let G be a graph with p vertices, and let L be its Laplacian matrix. Goemans and 
Williamson [Tl] introduced the following semidefinite programming-based upper bound on 
maxcut(G): 



maxc-uKCr) < QW{G) := max <j ^trace(LX) 



X h 0, Xu = l(l<i<p) 



(5) 



It was shown in p] that 0.878GW(G) < maxcut(G) < GW{G) holds for all graphs G. 
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The associated dual semidefinite program takes the form: 



gW(G) = min \ V 
wet 1 



w. 



Diag(w) - -L y 



(6) 



where Diag is the operator that maps a p-vector to a p x p diagonal matrix in the obvious 
way. 



3.1 The Goemans- Williamson bound for G n 

Using the technique of symmetry reduction for semidefinite programming (see e.g. [10J), 
one can simplify the dual problem ([6]) for the graphs G n defined in Section by using the 
dihedral automorphism group of G n . We state the final expression as the following lemma. 

Lemma 4. Let n > be an odd integer and d = \n/2\ . One has 



where 



min i jn^i/j Diag (^y - ^val^j + A (m) >z (0 < m < d) j 



vol; 



(m) 



A 



iyi 

1 

4 



l) + 2(i-l)(d-i), 2<i<d, 



i-l 



- 2ivrnk -J — 1 



n— 1 

+ E ■ 

fc=l k=n— j+1 



-27rmfc V-— T 



, 2<i<j<d, 



(7) 



For the proof, we recall that the Kronecker product A (g> B of matrices A G M pX9 and 
B € M rx,s is defined as the pr x qs matrix composed of pq blocks of size r x s, with block 
ij given by a^B where 1 < i < p and 1 < j < q. 

Proof. We first label the vertices G n as follows. Consider the cycle C n with vertices num- 
bered {0, 1, . . . , n — 1} in the usual way. The vertices of G n that correspond to chords con- 
necting points at cyclic distance i are now given consecutive labels (0, i), (1, i + 1), . . . (n — 
1, i — 1). Thus the adjacency matrix of G n is partitioned into a block structure, where each 
row of blocks is indexed by a cyclic distance i £ {2, . . . , d}, and each block has size n x n. 

Moreover, block € {2, . . . ,d}, i < j) is given by the n x n circulant matrix with 

first row 

[0 n _ i _ J - +1 Oj_j], 
where and denote the all-ones and all-zeroes vectors in R fc , respectively. 
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The eigenvalues of this block are 



e s + e s (0<m<n-l); (8) 

fc=l k=n— j+1 



SCO C 



•g. US]. 



Now let an optimal solution w of the semidefinite program (|SJ) be given for G = G n . If we 
project the matrix 

Diag(w) + -L 

onto the centralizer ring of Aut(G n ), then we again obtain an optimal solution. Indeed, this 
projection simply averages the components of w over the d— 1 orbits of Aut(G n ). Moreover, 
the projection is also a symmetric positive semidefinite matrix, since any projection of a 
Hermitian positive semidefinite matrix onto a matrix *-algebra is again positive semidefinite 
(see e.g. [II]). 



Denoting the average of the w components in orbit i by yi, we obtain an optimal solution 
of the form 

d 

QW(G n ) = min nS^yi 



i=2 

such that 

d 1 

Vi far-ieHi) ® In - y 0, (9) 

where e» denotes the i-th standard unit vector in and I n denotes the identity matrix 

of order n. 

Let Q denote the (unitary) discrete Fourier transform matrix of order n. Condition @ is 
equivalent to 

(I n (g) Q) fa Vi ( e i-l e f-l) ® !n - (In ® Q)* >= 0. (10) 

Since the unitary transform involving Q diagonalizes any circulant matrix (see e.g. [13J), 
the matrix (I n ® Q)L(I n x Q)* becomes a block matrix where each nx n block is diagonal, 
with diagonal entries of block given by the eigenvalues in (|8|). 



Finally, the rows and columns of the left hand side of (|10p may now be re-ordered to form a 
block diagonal matrix with n x n diagonal blocks given by the right hand side of (|7|) (only 
d + 1 of these blocks are distinct). This completes the proof. □ 

A few remarks on the semidefinite programming reformulation in Lemma HJ 

• The constraints involve Hermitian (complex) linear matrix inequalities, as opposed to 
the real symmetric linear matrix inequalities in ([6]). 
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• The reduced problem has d + 1 linear matrix inequalities involving (d — 1) x (d — 
1) matrices. By comparison, the original problem had one linear matrix inequality 
involving ((") — n) x ((") — n) matrices. As a result, the reformulation of QW(G n ) 
may be solved for much larger values of n than the original formulation ([6]) (see next 
section). 

• Although we have only done the symmetry reduction of problem © for G n with n 
odd, the case for even n is similar, but omitted, since we will not use it later. 

• Any feasible point y E of the reduced problem in Lemma 2] provides a certificate 
of an upper bound on QW(G n ), and consequently a certificate of a lower bound on 
u 2 (K n ), since u 2 (K n ) > - GW(G n ). 

4 Numerical computations: proof of Theorem [I] 

Theorem Q] (A) follows by an exact computation of the related maxcut problem of G n 
for certain values of n, while Theorem [1] (B) follows by a calculation of QW(G%§§) and a 
standard counting argument. 

4.1 Proof of (A) 

First we observe that if n < 5 then Z(n) = 0, and the assertion v<i{Kri) = Z(n) is easily 
verified. 

We computed the exact value maxcut(G n ) for n = 5,7,9,11,13,15,17,20, and 24, using 
the solver BiqMac [31], available from htt p : //biqmac . uni-klu . ac . at/| Computation was 
done on a quad-core 2.0 GHz Intel PC with 10 GB of RAM memory, running Linux. We 
used a cut-off time of 60 hours for the computation for each value of n. As a consequence, 
the BiqMac solver failed to terminate successfully in a few cases, namely n = 19, 21, 22, and 
23. 

The results are presented in the second column of Table |4*TT1 The exact value of V2{K n ) 
(fourth column) follows from the second and third columns (using Lemma [3]) . The fifth 
column is given for reference, to verify that viiKn) = Z{ n ) f° r an these values of n. Thus 
(A) follows for n = 5,7,9,11,13,15,17,20, and 24. The last two columns show the 
CPU time required, and the number of nodes evaluated in the branch an bound tree by 
the solver BiqMac. Finally, an elementary, well-known counting argument shows that if 
^(-K^m+i) = Z(2m + 1) for some positive integer m, then ^2(^201+2) = Z(2m + 2). This 
proves (A) for the remaining cases n = 6,8, 10, 12, 14, 16, and 18. 
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Table 1: The second column gives the exact values of maxcut(G ri ) that we computed. 
The fourth column gives the corresponding exact values of V2(K n ) (using that v^iKn) = 
\E n \ — maxcut(G n )). For all these values of n, the conjecture V2{Kn) = Z(n) is verified. 



4.2 Proof of (B) 

The first ingredient in the proof of (B) is a lower bound for ^2(^899)- We obtained this bound 
via the approximate calculation of C/W^Gsgg), which we achieved by using the semidefi- 
nite programming reformulation in Lemma [H Computation was done on a Dell Precision 
T7500 workstation with 92GB of RAM memory, using the semidefinite programming solver 
SDPT3 [351 [37] under Matlab 7 together with the Matlab package YALMIP [26] • The total 
running time was 12, 602 seconds. SDPT3 was chosen since it can deal with Hermitian 
matrix variables. We obtained ^>V(G 899 ) < 1.76537474 x 10 10 . Using Lemma[3]and ©, it 
follows immediately that 

^2(^899) > 9, 381, 181, 976. (11) 

The second ingredient to prove (B) is to establish a lower bound on the asymptotic ratio 
\\m n ^ too V2{K n ) / Z{n) that can be guaranteed from a lower bound on i>2{Km) for some 
m > 3. 

Claim 5. For any integer m > 3, 

11111 VI \ ~ —( TV OW 9T M K m)- 

n->oo Z(n) m{m — l){m — 2){m — 6) 

Proof. Let m, n be integers with 3 < m < n. Consider a 2-page drawing D of K n with 
V2{K n ) edge crossings. Let Q denote the set of subgraphs of K n that are isomorphic to K m , 
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i.e. \Q\ = Any two disjoint edges in K n occur in of the graphs in Q. Thus, every 

crossing in D appears in the induced drawings of ("Z4) graphs in Q. Consequently, 

(K s. V2(K m )Q = V2(K m )n(n-l)(n-2)(n-3) 
2[ n) ~ (™I 4 4 ) ro(m-l)(m-2)(m-3) ' 

The claim follows immediately from this inequality and the definition of Z{n). □ 

It only remains to observe that (B) is an immediate consequence of and Claim [SJ 



5 A quadratic programming lower bound for v^iKm^) 

Throughout this section, assume that m is fixed, and consider 2-page drawings of K m>n , 
where n is any positive integer. Thus, all vertices lie on the x-axis, and each edge is contained 
either in the upper or in the lower half-plane. We assume, without any loss of generality, 
that the m degree-n blue vertices bi, 62, . . . , b m appear on the x-axis in this order, from left 
to right. The n degree- m vertices are red. The star of a red vertex r (which we shall denote 
star(r)) is the subgraph induced by r and its incident edges. Thus, for every red vertex r, 
star(r) is isomorphic to K m ^\. 

5.1 The type of a red vertex 

In our quest for lower bounding the number of crossings in any 2-page drawing T> of K mjn , 
the strategy is to consider any two red vertices r, r' , and find a lower bound for the number 
x v (star(r), star(r')) of crossings in T> that involve one edge in star(r) and one edge in 
star(r'). The bound we establish is in terms of the types of r and r'. The type (formally 
defined shortly) of a red vertex is determined by its position relative to the blue vertices, 
and by which edges incident with it lie on each half-plane. 

We start by noting that we may focus our interest in drawings in which no red vertex lies to 
the left of b\. Indeed, if the leftmost red vertex lies to the left of 61 (and so it is the leftmost 
vertex overall), it is easy to see that it may be moved so that it becomes the rightmost 
(overall) vertex, without increasing the number of crossings. By repeating this procedure 
we get a drawing with the same number of crossings, and with no red vertex to the left of b\. 
Thus there is no loss of generality in dealing only with drawings that satisfy this property, 
and it follows that each red vertex r has a position p{r) relative to the blue points: p(r) is 
the largest j G {1, 2, . . . , m} such that r is to the right of bj. 

Also, to each red vertex r we can naturally assign a partition {U(r), L(r)} of {1,2,..., m}, 
the distribution of r, defined by the rule that j G {1,2,... ,m} is in U(r) (respectively, 
L(r)) if the edge rbj lies in the upper (respectively, lower) half-plane. We call the triple 
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(p(r),U(r), L(r)) the type of r, and denote it by type(r). Since p(r) can be any integer 
in {1,2,..., m}, and U (r) any subset of {1,2,..., m} (and L(r) = {1,2,..., m} \ U(r) is 
determined by U(r)), it follows that there are m2 m possible types for a red vertex. We use 
Types(m) to denote the collection of all m2 m possible types. 

5.2 Guaranteeing crossings between red stars using types 

The motivation for introducing the concept of type is that knowing the types of two red 
vertices r and r' in a drawing T> of K m ^ n yields a lower bound on x P (star (r), star (r')). We 
illustrate this with an example. Suppose that m = 5, and that type(r) = (2, {1, 2, 3, 5}, {4}) 
and type(r') = (3, {1, 3, 4, 5}, {2}). The situation is thus as illustrated in Figure HI 




x 



Figure 4: The types of the red vertices r and r' are (2, {1, 2, 3, 5}, {4}) and type(r') = 
(3, {1, 3, 4, 5}, {2}), respectively. Thus, r is in position 2 (that is, between 62 and 63), and 
the edges joining r to 61, 62, 63 and 65 are in the upper half-plane and the edge joining r to 
64 is in the lower half-plane. Both crossings in this drawing can be easily predicted from 
type(r) and type(r'). 

Both crossings between star(r) and star(r') in this example are easily detected from type(r) 
and type(r'). Indeed, since b%,r, r', 65 occur in this order from left to right (this follows 
since r and r' are in positions 2 and 3, respectively), and b\r' and rb§ are both on the upper 
half-plane (this follows since 1 € U(r') and 5 6 U(r)), it follows that b±r' and rb§ must 
cross. We remark that the key pieces of information are that (i) the endpoints b\,r, r' , 65 of 
b\ r' and rb^ alternate on the x-axis (that is, they are all distinct and occur in the x-axis so 
that the ends of one edge are in first and third place and the ends of the other edge are in 
second and fourth place); and (ii) both edges are drawn on the same half-plane. 

Using this simple criterion (if two edges are on the same half-plane and their endpoints 
alternate, then they must cross each other), given two red points r,r' in a drawing T> of 
K-m,n, it is easy to derive a lower bound for x B (star(r), star(r')) in terms of type(r) and 
type(r'). This bound (Proposition [6] below) is given in terms of a quantity we now proceed 
to define. 
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First, for a = (p, U, L) and r = (p 1 , U', V) G Types(m), we let 

: | (*>.?') | ({i G U and j G U') or (i G L and j G L')^ and 

(^(i < j < p) or (j < p and p' < i) or (i < j and p < i) or (p < j < i < p) ^ | 



<t,t := 



and 



'[cr,r], if p<p', 

[r,cr], if p>p', 

,min{[(T,r], [r,cr]}, if p = p'. 



The nonnegative integers Qo-t- can be naturally regarded as the entries of a m2 m x m2 m - 
matrix Q indexed (both by rows and columns) by the elements of Types(m). It is easy to 
check that the matrix Q is symmetric, and its entries provide the lower bounds we have 
been aiming for. 

Proposition 6. Let a,r G Types(m), and suppose that r a ,r T are red points in a drawing 
V of K m ^ n , such that type(r CT ) = a and type(r r ) = r. Then 

x c (star(r (T ),star(r r )) > Q aT . 

Proof. Suppose first that r a occurs to the left of r T . It is easy to verify that if i,j are 
integers such that either (i) i < j < p; (ii) j < p and p' < i; or (hi) i < j and p' < i; or 
(iv) p < j < i < p' , then the endpoints of rbi and r'bj alternate. Therefore, if either i G U 
and j G U', or i G L and j G V , then rbi and r'bj cross each other. Therefore there is an 
injection from the set of all pairs (i,j) of integers that satisfy the condition in the definition 
of [a, r], to the set of crossings that involve an edge in star(rv) and an edge in star(r T ); that 
is, x x ,(star(r CT ), star(r T )) > [a,r]. 

Similarly, if r a occurs to the right of r T , then x c (star(r (7 ), star(r T )) > [r, a]. 

Now if p < p' (respectively, p > p'), then r a necessarily occurs to the left (respectively, 
to the right) of r T , and so it follows that x E ,(star(r (T ), star(ry)) > [a, r] = Q aT (respec- 
tively, > [r, a] = Qctt)-, as required. Finally, If p = p', then r a can be either to the 
right or to the left of r T . In the first case, x B (star(rv), star(r T )) > [a, r], while in the 
second case x E ,(star(r (7 ), star(r T )) > [t,o~]. Thus, in this case, x c (star(r cr ), star(r T )) > 
min{[<7, r], [r, a]} = Q aT , as required. □ 



5.3 The quadratic program 

Consider now any fixed 2-page drawing V of K m ^ n . For each type a G Types(m), let no- 
denote the number of red vertices whose type in V is a, let p a := n a /n, and let p be the 
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vector (Po-)crgTypes(m)- It follows immediately from Proposition [6] that the number U2(T>) of 
crossings in T> satisfies 

(P) > 2 £ o 

(TT f7(T I .x I 

o-,reTypes(m) (jGT yP cs(m) 

CT,rGTypes(m) <T£Types(m) 

= y P QP - g 2^ 0<t«tP«t 

(XGTypes(m) 

n 2 T n 
> ttP Qp - o max , s ^ 

Z Z cr£Typcs(m) 

n 2 t.„ m(m — l)n 
— ~2" P ^ P 4" 

where for the last inequality we use that S . gTypcs ( m )p (J = 1 and that Qo-o- = [a, a] < (™) . 

The derived inequality holds for every 2-page drawing T> of K mjn , and so in particular for 
a crossing-minimal drawing. Thus, if we let 



A= |x= (xi,x 2 ,...,x m2 ™) T £l mr 
denote the standard simplex in R m2m , then we obtain 



y^gj = 1, Xj > oj 



n 2 f . r \ m(m — l)n 
^ 2 (^ m ,n) > -TT mm x T Qx - V ; . (12) 
z yxeA y 4 

We may therefore obtain a lower bound on v<2,{K m ,n) f° r some fixed m (we will be particularly 
interested in the case m = 7), by solving the standard quadratic programming problem 

lb(m) = min x T Qx. (13) 
xeA 

The standard quadratic programming problem is NP-hard in general, and we will only 
compute a lower bound on the minimum via semidefinite programming, as explained in the 
next section. 
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6 A semidefinite programming lower bound on iS2(K m ,n) 

The usual semidefinite programming relaxation of problem (|13p takes the form 

lb(m) > min{trace(QX) | trace(JX) = 1, X y 0, X > 0} 
= max{t | Q - tJ = Si + S 2 , Si y 0, S 2 > 0} 

:= SDP 6oW (m), (14) 

where J is the all-ones matrix, and X > means that X is entrywise nonnegative. We 
observe that the first equality is due to the duality theory of semidefinite programming. 

Due to the special structure of Q, we may again use symmetry reduction to reduce the 
size of these problems. To this end, for odd m, we may order the rows and columns of 
Q to obtain a block matrix consisting of circulant blocks of order 2m. (Thus there are 
2 m_1 rows/columns of blocks). The ordering of rows works as follows: we first define a 
group action on the set Types(m). For ease of notation we now represent the elements of 
Types(m) as (p, U), with p G {0, . . . , m — 1} and U C {0, . . . , m — 1}, i.e. we now number 
the m vertices from to m — 1, and omit the set L (which is redundant in the description 
since it is the complement of U). 

The group in question is generated by the following two elements, a 'flip': 

<7i : (p,U) ^ (p,{0,...,m-l}\U), 

and a 'cyclic shift': 

92 '■ (p, U) t— > (p + 1 mod m, {u + 1 mod m \ u € U}). 

Note that gi and 52 commute and therefore generate an Abelian group of order 2m. If m 
is odd, then g := g% o g 2 generates the entire group, i.e. in this case we obtain the cyclic 
group of order 2m. Indeed, the order of g equals the least common multiple of the orders 
of g\ and g 2 , namely 2m if m is odd. 

Also note that 

Qa,r = Qg,(a),g,(r) Vcr, T G Types(m), I G {1, 2}, 

i.e. the crossing number of a 2-page drawing does not change if we 'flip' the drawing along 
its spine, or, in the circular model, rotate the drawing. 

Finally, we group together the 2m elements of Types(m) that belong to a given orbit of the 
group, to obtain 2m x 2m circulant blocks. In what follows, we denote the first row of the 
2m x 2m circulant block (i, j) by q^'^ G l? m . 

Lemma 7. For odd m, the semidefinite programming bound |i^[ ) may be reformulated as 

SBP bound (m) = maxt 
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subject to 

qf j) - t - xf j) > 0, < k < 2m - 1, 1 < i,j < 2 m ~ 1 , 

2m- 1 



Xf> = xf j) + Y, xf 3) e-^ uk ' m , l<i<j< 2 m -\ < t < 2m - 1, 

k=l 

x (t) = ( X (t)y y < i < 2m - 1, 

4 M) -4t?+l-fc = °' l<fc<m-l, l<i<2 m -\ 
x (i,j) € R 2m^ l<iJ<2 m -\ 

Proof. The proof is similar to that of Lemma H] and is therefore omitted. □ 
A few remarks on the semidefinite programming reformulation in Lemma [7] 

• As in Lemmadl the constraints involve Hermitian (complex) linear matrix inequalities. 

• The reduced problem has 2m linear matrix inequalities involving (2 m_1 ) x (2 m_1 ) ma- 
trices. By comparison, the original problem had one linear matrix inequality involving 
a (m2 m ) x (m2 m ) nonnegative matrix. As a result, the reformulation in Lemma 
may be solved for larger values of m than the original formulation (see next section). 

• Similarly to Lemma [H every feasible point x^'^ € M? m (1 < i,j < 2 m_1 ) yields a 
certificate of lower bound on SDPft oun( ^(m), and consequently a certificate of a lower 
bound on u 2 (K mtn ), by (H2]). 



7 Numerical computations: proof of Theorem [2] 

Using the reformulation in Lemma we showed numerically that SDPb otm d(7) = |. Com- 
putation was done on a Dell Precision T7500 workstation with 92GB of RAM memory, 
using the semidefinite programming solver SDPT3 |35|. [37] under Matlab 7 together with 
the Matlab package YALMIP [26]. The running time was 23,774 seconds. SDPT3 was 
chosen since it can deal with Hermitian matrix variables. 

Using that SDP boW (7) = 9/2, it follows from CO]), CG|, and CHD that 

MK 7 ,n) > (9/4)n 2 - (21/2)n. (15) 

We recall that Z(7,n) = 9[n/2\ [(n - 1)/2J = (9/4)n 2 + O(n), and that v 2 {K 7 , n ) < Z(7,n) 
(since there are 2-page drawings of Kj tn with exactly Z(7,n) crossings). Using these obser- 
vations and (|15p , Theorem [2] follows for m = 7. 

Now an elementary counting argument shows that V2{K%, n ) > 8^2(^7^ )/6, and so using 
(|15p and simplifying we obtain ^2(^8, n) > 3re 2 — 14n. Since Z(8, n) = 3n 2 + 0(n), Theorem 
1 follows for m = 8. 
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8 Concluding remarks 



The Goemans- Williamson bound (Section [3]) empirically yields better lower bounds on 
V2{K n ) as n grows; see Figure [HI 



T 




I I I I I I I I I 

100 200 300 400 500 600 700 800 900 



Figure 5: The ratio j^f Gn) for n = 99, 199, 299, 399, 499, 599, 699, 799, and 899. 

Based on this empirical evidence, it seems reasonable to expect that the constant 0.9253 
would be improved if GW(G m ) were computed for larger values of m. Having said that, 
the figure also shows a trend of diminishing returns — by extrapolating the curve in the 
figure, it seems that it may not be possible to improve the constant to more than 0.929, 
say, through computation of GW(G m ), if m < 2,000. 

Another possibility to improve the constant is to compute i>2(K m ) for larger values of m 
than m = 24, by solving the maximum cut problem in Lemma El If, for example, one 
could verify in this way that ^2(^30) = ^(30), then this would yield the constant 0.9297, 
by Claim [5j 

Regarding the computational lower bound on V2{K m , n ): It is interesting to note that the 
SDP bound SDP bound(rn) provided a tight asymptotic bound on z^O^m,™) for m = 3, 5 and 
7. A similar SDP bound used in [24] and [23] did not provide a tight asymptotic bound on 
the usual crossing number cT(K mtn ), not even for m = 5. Our results therefore suggest that 

one may be able to prove computationally that limn^oo V2 ^m^n) = ^ ^ or (fi xec odd values 
of m > 9. Having said that, for m = 9, the resulting semidennite program was too large 
for us to compute SDPf, onra d(9). This problem therefore provides a good future challenge to 
the computational SDP community. 
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